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ABSTRACT 

We present a fully parallel chemo-dynamical Tree/SPH code, GEAR, which allows to perform high resolution simulations with detailed 
chemical diagnostics. Starting from the public version of Gadget-2, we included the complex treatment of the baryon physics: gas 
cooling, star formation law, chemical evolution and supernovae feedback. We qualified the performances of GEAR with the case of dSph 
galaxies. GEAR conserves the total energy budget of the systems to better than 5% over 14 Gyr and proved excellent convergence of the 
results with numerical resolution. We showed that models of dSphs in a static Euclidean space, where the expansion of the universe is 
neglected are valid. In addition, we tackled some of the existing open questions in the field, like the stellar mass fraction of dSphs and 
its link with the predicted dark matter halo mass function, the effect of the supernova feedback, the spatial distribution of the stellar 
populations, and the origin of the diversity in star formation histories and chemical abundance patterns. Strong supernovae driven 
winds seem incompatible with the observed metallicities and luminosities. Despite the fact that newly formed stars are preferentially 
found in the galaxy central parts, turbulent motions in the gas can quickly erase any metallicity gradient. The variety in dSph properties 
result from a range of total masses as well as from a dispersion in central densities. The latter is also seen in the haloes emerging from 
a ACDM cosmogony. 
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1. Introduction 

The first motivation of the present work was the clear need for 
an extension of ^Revaz et al. ( 2009 1, which sought to identify 



the dominant physical processes at the origin of the dynamical 
and chemical properties of dwarf spheroidal galaxies (dSphs). 
While allowing a critical step forward with a large set of simula- 
tions, the excellence of the treeAsph code used in this work was 



somewhat lessen by its serial numerical structure (Serna et al. 
|1996[[Aiimi et al. 2003 ), hampering the possibility of high reso- 
lution simulations due to large computing time. 

The analysis of other Tree/SPH codes, that include chem- 
ical evolution schemes and were applied to the formation and 
evolution of dwarf galaxies, demonstrates that the specific astro- 
physical aims of each study influence if not drive the numerical 
implementations. Using colour-magnitude diagrams as observa- 
tional probes, Carraro et al. ( 200 1| modeled their gas particles 
as closed box models and used the instantaneous recycling ap- 
proximation. Marcolini et al. (2006, 2008) focused on the re- 
sponse of the ISM to the supernovae explosion and introduced 
an arb itrary num ber of star formation episo des. [Kawata et al.| 
( |200 6^ following Kawata & Gibson (2003'), combined kinetic 
and thermal feedbacks, which could generate galactic winds in 
elliptical galaxies. |Okamoto et al.| ( [2010| l investigated the effect 
of different recipes of kinetic feedback driving galactic winds 
on the global properties of the Milky Way satellite model galax- 
ies. [Read et aL] (|2006') looked for a mass boundary distinguish- 
ing between essentially pure dark haloes and small stellar sys- 
tems. They chose a unique global yield reproducing the metal- 
licity of the inter-galactic medium at z=3. |Stinson et al.] ( |2007| l 



and [Governato et al.| ( [20T0| l applied the 'S tinson et al.' ('2006Vs 
blast wave recipe for the feedback. Sawala et al. " ^2010 ) used 
the multiphase scheme for the interstellar medium developed by 
Scannapieco et al. (2005 2006 1, and constrained the fraction of 
metals given to the cold gas phase, later defining the final metal- 
licity for a given stellar mass, on the metallicity-relation of dwarf 
galaxies. Valcke et al.| ( |2008) and later Schroyen et al. ( 201 l| l 
took into account stellar winds and supernovae explosions but 
only consider the global metallicity Z. Ricotti & Gnedin ( 2005| l 
stopped their simulations very early on (z~8) but included con- 
tinuum radiative transfer to compare the global properties of 
dSphs with a simulated sample of galaxies formed before reion- 
ization. 

We aimed at a new N-body code, hereafter GEAR, which ful- 
fills a number of requirements among which piimaiily the ca- 
pability to offer high spatial resolution together with detailed 
chemical diagnostics, and to follow the galaxy evolution over a 
full Hubble time, either in isolation or in a cosmological context. 
Because our approach to constrain the scenarios of galaxy evolu- 
tion involves the confrontation with observed stellar metallicity 
distributions and stellar abundance ratios, we seriously evaluated 
the impact of the supernovae feedback on these quantities, and 
further probed all ingredients in the parametrization of the star 
formation. While it would be illusive to expect tracing in detail 
the star formation processes at the scale of molecular clouds, still 
we tried to release as much as possible ad hoc assumptions, and 
based our modeling on physical grounds. Great care was devoted 
to trace all sources and sinks of energy and guaranty an accurate 
budget of gains and losses. We also controUed the convergence 
of the results. 
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GEAR is a flexible fully parallel chemo-dynamical Tree/SPH 
code, which we applied to tackle some of the questions left open 
after [Revaz et al. (j2009 ), such as the gas motions and the final 
radial distribution of the stellar population in dwarf spheroidal 
galaxies. It also allowed to expand on a number of new subjects, 
in particular the question of the relation between galaxies in iso- 
lation and their parent haloes in cosmological simulations, and 
clarify the final baryonic fraction of galaxies as a function of 
their halo masses. 

The paper is organized as follows: the code and the imple- 
mentation of physical processes are described in Sect. |2] The 
initial conditions are detailed in Sect. |3] The robustness of the 
code is discussed in Sect. |4] The role and impact of the main 
parameters are evaluated in Sect.|5] while Sect.|6]focuses on the 
initial mass and central gas density. Sect. |7] describes the differ- 
ent regimes of star formation for low mass galaxies and makes 
the first detailed comparison with the observations. Sect. |8] ad- 
dresses the question of the spatial distribution of the stellar pop- 
ulation as a function of their age and metallicity. Sect. |9]investi- 
gates the link between stellar masses in galaxies and their parent 
dark matter haloes. Sect. [Tolsummarizes our results. 



2. GEAR 



We supplemented the public version of Gadget-2 (Springel 



2005 j ) with the complex treatment of the baryon physics: gas 
cooling, star formation law, chemical evolution and supernovae 
feedback. 



2.1. The gas cooling 

The interstellar medium (ISM) is modeled as an ideal, inviscid 
gas with an adiabatic index y = 5/3. The continuity equation 
is inserted into the first law of thermodynamics in order to fol- 
low the evolution of the gas specific internal energy. The varia- 
tion of the internal energy depends on the mechanical forces, the 
artificial viscosity and the radiative cooling expressed through 
a cooling function A(p, T). Following the integration scheme 
of Gadget-2, we use the entropy function A (see Springel & 



[Her nquist (2002) for the exact definition) instead of the internal 
energy as independent thermodynamic variable. The variation 
of A is simplified as a sum of two terms corresponding to the 
contributions of the artificial viscosity ( Springel||2005 1 and the 
cooling: 



d , d , 
at at 



where: 



dt 



cool 



dt 



-{y-l)p-^K(p, T). 



(1) 



(2) 



For temperatures hotter than 10'* K, the cooling function is 
calculated following the me tallicitv dependent prescription of 
Suflierland & Dopita| ( [l993 ). Below lO'^K, the cooling of H2 
and HD molecules are taken into account as well as the atoms of 
oxygen, carbon, silicon and iron ( |Maio et aL]|2007| l, see Revaz | 
|etal.| ( |2009| ) for more details. 

In dense (p > O.Olme/cm^) and warm-hot regions {T > 
10^ K), the cooling time may be much shorter than the dynamical 
time requiring extremely short time steps. In such a case, we 



use the isochoric approximation (see for example, Springel et al. 
( |2001| l), assuming a constant density p. The integration of Eq. []_ 



is performed using adaptive time steps set to a fraction of 2% of 
the cooling time defined as : 



^cool 



(3) 



cool 



We do not take into account the heating of the gas by a cosmic 
background UV field. Its effect was investigated in a number of 
former studies (e.g, Sa wala et al.| (2010); Okamoto et al. (2010)), 
because it could contribute to evaporate the gas in small haloes 
(vmax < 12km/s). 

However, the ISM is known to be strongly non homogeneous 
down to very small scales, reaching densities ~ lO^mn/cm-* 
(e.g., Omont 2007; Be uther et al.|2007 and references therein), 
for which the cooling time is very short. As long as those over- 
densities are not properly resolved, including an additional heat- 
ing source, as the UV background is not a real gain, because the 
cooling will be underestimated due to the lack of resolution. This 
is without taking into account that the gas may be self-shielded 
strongly complicating the effect of the UV field. In this picture, 
improving the cooling by following precisely all cooling agents 
involved in is also unprofitable. As no heating process is taken 
into account here, the integration of Eq. [T] does not suffer from 
numerical instabilities and an implicit integration scheme is not 
required. 



2.2. Star formation iaw 

The complex and poorly understood formation processes is im- 
plemented using the phenomenological prescription proposed by 
[Katz (1992 ) and |Katz et aL] fl996). 

A gas particle becomes eligible to star formation when it ful- 
fills the following physical conditions: /) the particle is collaps- 
ing (its velocity divergence is negative), //) its density is higher 
than a threshold psfr, Hi) its temperature is lower than a threshold 

7"sfr- 

In a time interval At, an eligible gas particle of ma ss nij, has 
a probability pi, to form a stellar particle of mass m*. ( Springel 
& Hemquist (2003)): 



nit, 



P* = 



1 - exp|-— Af 



(4) 



where c* is the star formation parameter and t„ the local free fall 
time. This ensures a star formation law independent of both the 
time step Af and m^, : 



dp* 
dt 



to 



(5) 



in which pg is the gas density. 

Each gas particle can form a number of stellar particles, this 
number is set by the parameter A^*. Each new stellar particle 
represents an ensemble of stars sampling an initial mass func- 
tion (IMF) with different slopes, a, in four stellar mass ranges 
(Kroupa 2001), taking into account the systematic effects due to 
unresolved binaries : 



-0.3, if m< 0.08 Mq 

-1.8, if 0.08 Mo <m< 0.50 Mq 
-2.7, if 0.50 Mo <m< 1.00 Mo 
-2.3, if m> 1.00 Mo 



(6) 



The minimal and maximal stellar masses, that are considered 
in the IMF, are 0.05 and 50 Mo, respectively. The stellar parti- 
cles receive the positions, velocities, and chemical abundances 
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of their parent gas particles. These gas and stars are dynami- 
cally decoupled, therefore their evolution in the phase space are 
quickly diverging. 

In GEAR, the stellar particles are numerically considered as 
an independent class of particles with their own structure. This 
is necessary in order to optimize the memory needed to store 
their chemical properties. GEAR keeps the possibility to order all 
types of particles following the Peano-Hilbert curve similarly 
to Gadget-2. This allows to take advantage of the fast cache- 
memory available on modern processors. 

We set Tsfr to 3 x lO^K. The exact value of T^tr is not cru- 
cial. It is chosen here to avoid the formation of stars in hot gas. 
The other parameters, c^,, psfr, and the maximal number of stellar 
particles that may be created from one gas particle A^* = Wg/m* 
are discussed in detail in Section|5] 

2.3. Chemical evolution 

The galactic evolution is influenced by two main features of the 
stellar evolution: on the one hand the nucleosynthesis, on the 
other hand, the stellar feedback (winds or supernova explosions). 
In the following, we only consider the effect of supernovae. 
Indeed, stellar winds from intermediate mass stars inject little 
power to the ISM as compared to the explosion of supernovae 
( [Leitherer et al.|[r992[ ). Moreover, the evolution of the galactic 
systems that we are considering is already well constrained by 
the abundances of or-elements (magnesium) and iron and their 
ratios. Both types of elements are produced by Type la (SNela) 
and of Type II (SNell) supernovae. 

For a given stellar particle, the number of stars ending their 
lifetime is computed at each dynamical time step. These stellar 
lifetime depends on the metaUicity at the creation of the parti- 
cle and are taken from Kodama & Arimoto ( 1997 private com- 



efFect of varying A^ngb is discussed in Section 5.6 In order to 



munication). For an accurate resolution, we impose a maximum 
time step of 0.2 Myr. This corresponds to less than 10% of the 
smallest lifetime of the exploding stars in our simulation (about 
3 Myr). The stellar mass dependent yields of SNell for stars be- 
tween 8 and 50 Mq are taken from Tsujimoto et al.'(T995). We 
adopt the model of Kobayashi et al. (2000| for the explosion of 
SNela. The progenitors of SNela have main-sequence masses 
between 3 and 8 Mq, and evolve into C-hO white dwarfs (WDs). 
These white dwarfs can form two different types of binary sys- 
tems, either with main sequence stars or with red giants (see 
|Revaz et al.| ( |2009| ) for more the details on the adopted parame- 
ters). The nucleosynthesis products of SNela are taken from the 



model W7 of Iwamoto et al. ( 1999 1. Contrary to Kobayashi et al. 
( |2000| l we do not prevent the SNela explosion at [Fe/H] below 
-1. We note that the nucleosynthesis tables and the IMF can be 
changed easily in GEAR . 

The stellar chemical ejecta are distributed among the nearest 
neighbors following the SPH scheme (e.g., Wiersma et al.'2009j 
Revaz et al.|[2009 ). First, the gas density p,- and corresponding 
smoothing length hi are determined at the location of each stellar 
particle Both quantities are computed similarly to the case of 
the gas particles (Springel 2005). The mass of metals ejected 
from particle / and attributed to particle j is computed through 
the weights Wy-: 



improve the conservation of energy during the process of mass 
ejection, the velocities and the entropy of the particles are mod- 
ified as described in Appendix [A] 

The chemical abundances are calculated with respect to the 
solar abundances of |Anders & Grevesse| ( |1989] l. The V-band lu- 
minosities are derived following the stellar population synthe- 
sis model of Vazdekis et al. (1996i computed for the revised 
(Kroupa 2001|'s IMF. Where necessary, the luminosities are 
inter- and extra- polated in age an metallicity using a bi variate 
spline. 

2.4. Feedback 

Each explosion of supernova injects an energy esN^SN into the 
ISM, with /isN - 10^' ergs and esN the feedback efliciency fac- 
tor, which is further discussed in Section [5T| 

The supernova explosion feedback is definitively a complex 
mechanism, which led to a large number of reports discussing 
both its numerical implementation and its physical operation 
modes. One can broadly define three different types of meth- 
ods discussed in the literature: /) The so-called thermal feedback 
in which the energy is released through heating of the gas : ei- 
ther in an homogeneous ISM (['Gerritsen|1997[|Mori et al.|1997 



Thacker & Couchman 2000 ; Sommer-Larsen et al. 2003 ; Brook 
et al. 2004, Stinson et al. 2006) o r in a m ultiphase ISM (Yepes 
l et al.l|1997'['|Hultman & Pharasyn||T999| |Marri & White||2003 ' 



'Scannapieco et al. 120061). ii) A kinetic feedback in which energy 
is mechanically released (Navarro & W hite 1993, Springel *&] 
Hernquist 2003] |DaUa Vecchia & Schaye|20 08). 

We were very careful in ensuring an accurate budget of gains 
and losses of energy, even in presence of a strong feedback. We 
also checked that the energy from the supernovae, as predicted 
by stellar evolution, was fully injected in the galactic system. 
This particular point is not always granted, especially when us- 



ing a stochastic implementation of the kinetic feedback ( Springel 
| & Hernq uist 2003 , |Dalla Vecchia & Sch^[2008l l. In practice, 
at each time step [f, f -H Af], we calculate the amount of energy do 
be released during the explosions of SNe. This quantity £'sn(0 
is the sum of the contributions of SNela and SNell: 



EsNiO = ^ AE,- SN = ^ '«;,*,o [nai(0 £11 + n/,ia(0 ^la] ■. 



(8) 



mjW(rij,hi) 



(7) 



where W is the SPH kernel function, nij is the mass of the 
particle j and r,j the distance between the particles / and j. By 
default, the number of neighbors used is set to A^ngb = 50. The 



where m, * is the initial mass of the stellar particle /, and n,,n(f) 
and n,,ia(f) are the corresponding numbers of supernovae SNell 
and SNela per unit mass during At. The integral of £'sn(0 over 
time is then compared to the energy injected in the system during 
the feedback procedure. 

In the case of a kinetic feedback, the particles are pushed 
away, forming a wind with a velocity of the order of lOOkm/s 
( [Springel & Hernquist] [2003] |DaUa Vecchia & Sch aye"2008]l. 
This hardly conserves energy to a satisfactory level. Energy bud- 
get is improved if the wind particles are decoupled from the other 
gas particles or if the dynamical time steps are shorten before the 
advent of the wind. Unfortunately, this latter solution increases 
significantly the computation time. Moreover, the choice of a 
given wind speed limits the amount of energy deposited in each 
neighboring particle. Hence the total feedback energy derived 
from stellar evolution is not always fully transferred, especially 
for slow winds. Consequently, this technique implicitly imposes 
strong winds. 

For all these reasons, we chose a thermal feedback. At each 
time step, the energy Zssn is distributed in thermal form among 
the surrounding particles using the same weight Wij as for the 
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distribution of the chemical elements (Eq. |7|. When combined 
with radiative cooling, such a thermal feedback can be ineffi- 
cient. Indeed, the radiative cooling of the heated gas becomes 
very large and the injected energy is instantaneously lost ( |Katz| 
T992|. In GEAR, the cooling of the gas particle receiving the ther- 
mal feedback is switched off for a short adiabatic period of time, 
fad (a few Myr) ( GeiTi tsen|[T997l |Mori et al l[T997l [Thacker &| 
Couchm an 2000; So mmer-Larsen et al.|2003[|Brook et al.|2004[ 
Stinson et al. 20061. As both supernovae types release the same 
amount of energy, the switch is applied to gas particles receiving 
feedback from both SNell and SNela, without distinction. Thus, 
the feedback from SNIa may be more effective than the one used 
by other authors ( Stinson et al.,2006 , for example). The impact 
of 4d is discussed in Section 5.5 We carefully checked that our 
implementation of feedback combined with the cooling led to 
converging results with decreasing time step and increasing res- 
olution. 

We did not use any time step limiter to prevent possible nu- 
merical problems induced by a pre-shock timestep too long com- 
pared to the shock timescale. This time step limiter is important 
when large differences in temperature occur between the pre- 
and post- shock media ( [Saitoh &^l akino 2009; Merli n et al.| 
2010] ). It is not critical in our case, as the net heating of the ISM 
due to SNe explosions never exceeds a factor 20. 



3. Initial conditions 

3.1. Evolution of small dark matter haloes 

We ran a ACDM cosmological simulation in order to study the 
profiles of dark haloes with masses between 10^ and 10^ Mg, 



that are typical of dSphs ( [Walker et aLl|2007 
[20081. The volume of the simulation is 2^ h ' 



[Battagha et al. 
Mpc^ and con- 

tains 134'217'728 dark matter particles which results in a parti- 
cle mass of 4.6xlO^M0/h and a softening length of 150pc/h. We 
took the cosmological parameters from the concordance ACDM 
flat universe based on WMAP V data combined with the Baryon 
Acoustic Oscillations (BAO) in the distribution of galaxies and 
distances measurements from Type la supernovae (SN) observa- 
tions Pinshaw et al. 2009, Ko matsu et al.|2009| l : Q,„ = 0.279, 
Qa - 0.721 and ho = 0.7. We extracted the dark haloes at differ- 
ent redshifts using the HOP algorithm ( [Eisenstein & Hut|1998| l. 
The halo mass is defined as the mass inside raoo, the radius at 
which the density of matter is 200 times the critical density of the 
Universe; the halo density profiles have a NFW shape ( [Navarro[ 
et al. 1996 1997| l. We found 144 haloes with masses in the range 
from 10° to 10^ Mq. The analysis of these haloes showed two 
main features: 

(i) for 50% of the haloes, the density profiles in physical co- 
ordinates (as opposed to comoving coordinates) are in place at a 
redshift of 6, i.e., at about 1 Gyr (see Fig[T]i. Whilst the Universe 
was in average about 350 times denser at z = 6 than at z = (in 
physical coordinates), the density profiles of the dark haloes are 
similar. 

(ii) Varying masses by a factor 10, from 10^ to 10^ M©, scales 
the density profiles by a factor 3-4 only (see Fig|2]i. 

The consequences for the evolution of dSphs, are (i) the ex- 
pansion of the universe leads to the formation of stable sys- 
tems after z = 6. In turn, this justifies models of dSphs in a 
static Euclidean space, where the expansion of the universe is 
neglected. The physics of baryons that depends on the density 
in physical coordinates (for example the cooling of the gas) is 
correct. It would not be the case if the density in the haloes 
would increase with the mean density of the universe at higher 




Radius [kpc] 

Fig. 1. Evolution of the dark halo density profile for a halo mass 
of 6 X 10*^ Mq as a function of redshift. The dashed line indicates 
the limit of the resolution corresponding to the gravitational soft- 
ening. 




Ratiius [kpc] 



Fig. 2. The density profiles at z = of 142 dark matter haloes 
with masses between 10^ and 10^ Mq. The dashed line indicates 
the limit of the resolution corresponding to the gravitational soft- 
ening. The gray surface show the domain coved by the pseudo- 
isothermal profiles (Eq. |9]l adopted in our simulations, for the 
minimum and maximum values of central densities and radii 
('"max)- Only two haloes among the 144 extracted strongly de- 
viate from the NFW profile and are not seen in this plot. The 
reason is that at z = 0, they both experiment a strong merging. 



redshift. (ii) 50% of the haloes experience only minor mergers 
since z = 6. Their profiles are stable and can be modeled as iso- 
lated systems, (iii) Although the densities of haloes with mass 
between 10** and 10^ Mq are very similar, still they exhibit a 
small dispersion, a factor 3 to 4, which can help understanding 
the variety in the observed properties of the dSph galaxies. 

3.2. Isolated systems 

The previous conclusions warrant the simulation of dSphs as iso- 
lated systems, hence neglecting the expansion of the Universe. 
This approach is convenient in terms of CPU time, allowing to 
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run a large number of simulations to explore a wide range of 
parameters, as well as to reach very high resolutions. 

Our isolated systems are initially spherical and contain dark 
matter as well as 15% of baryons in form of gas. In a first ap- 
proximation, we assume the same profile for the gas and the dark 
matter. 

Instead of using an NFW profile, we preferred a cored one, 
which is supported by observations of normal low brightness and 
dwarf galaxies (Blais-Ouellette et ar||2001 de Blok & Bosma 



200?, 'Swaters et al. "2 003 [ [Gentile et al."20d4','200 5||Spekkens 
et al. 2005; de Blok 2 005[|de Bloket al . 2008; Spa no et al.|2008 



Walker & Peflarrubia|201 1 1, including the recent and high reso- 
lution observations of the THINGS survey ( |Oh et al.poTT) . Both 
gas and dark matter follow a pseudo-isothermal profile: 



Piir) = 



2' 



(9) 



where r is the radius, is the scale length of the mass distribu- 
tion, and pci the central mass density of the component /. 

All parameters are set to provide profiles that are compatible 
with the ACDM cosmological simulation discussed above, ex- 
cept in the central regions, where the profiles are flattened, see 
the gray area of Fig.|2] 

Both the initial dark matter and gas profiles are truncated at 
Tmax and the total initial mass is then defined as : 



a 



Tota] 

Tnternal 

PoteDtial 

Kinetic 

Cooling 





Time [Gyr] 



Fig. 3. Global energy budget for a dense and massive model of 
9.5 X IO'^Mq (p,,gas=0.071mH/cm\ = 0.05 and esN = 0.03) 
as a function of time. The total energy (in black) is the sum of 
all other lines, namely:the total potential energy (in red), the to- 
tal kinetic energy (in green), the gas plus stellar internal energy 
(in yellow), the radiative cooling energy (in blue), and the SNe 
feedback energy (in orange). The latter is negative because it 
corresponds to energy injected into the system. 



■ 47rpc,t 



arctan 



(10) 



Withpctot - Pchalo +Pc,gas- 

As the precise value of has only a limited influence on the 
evolution of the dSphs ( |Revaz et al.|2009 |l, we fixed it to 1 kpc. 
Pc.haio ™d Pcgas ^re related by the baryonic fraction. 

In the case of spherical systems and for an isotropic veloc- 
ity dispersion, the velocities of a component / can be derived 
using the the second moment of the Jeans equation (Binney & 
Tremaine] 1 987] |Hernquist]T993) : 



1 r 

Piir) J, 



dr'pi(r')dr'mr'). 



(11) 



For the halo, the velocity dispersion is directly taken from 
Eq.— 
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For the gas, the initial velocities are set to zero but its 
initial temperature is obtained by converting the kinetic energy 
of Eq. 1 1 into thermal energy : 



T(r) = /vu(r ■ 



(12) 



The parameter /vii set to 0.5 allows the gas to slowly flow 
towards the center of the system; , niu is the hydrogen mass, 
the mean atomic mass of the gas and ks the Boltzman constant. 



4. Robustness 

4.1. Energy conservation 

The first and fundamental requirement of numerical simulations 
is to carefully trace all sources (gain and loss) of energy in the 
systems : the total potential energy, the total kinetic energy, the 
gas internal energy, the radiative cooling energy and the SNe 
feedback energy. We also include in the gas internal energy the 



internal energy of gas particles converted into stars. The preci- 
sion of our energy budget is always better than 5%. As an ex- 
ample. Fig. [3] displays the global budget of energy of a mas- 
sive 9.5 X W^Me system, with = 0.05 and esN - 0.03, 
Pc,gas=0.071mH/cm-'. The total number of particles is 524'288 
and the softening length 25 pc. Such a high mass system, with a 
high central mass density, experiences star formation sustained 
at a high level (see Section |7]l, making it a difficult case from a 
numerical point of view. Nevertheless, the energy is conserved 
within 5%. This is remarkable with regards to the long integra- 
tion time (14 Gyr) compared to the dynamical time of the system 
(about 50Myr). Fig. [3] also reveals a clear balance between the 
injected energy by the SNe explosions and the radiative cool- 
ing. This means that although the cooling is switched off during 
an adiabatic period t.^^ (see Section 2.4 1, a large fraction of the 
feedback energy is still radiated away. 

4.2. Convergence 

We checked the convergence of our results with increasing res- 
olution, i.e., the number of particles was increased and the grav- 
itational softening eg was decreased accordingly. For this, we 
performed height sets of simulations on systems with two differ- 
ent masses, 3 and 9.5 x lO^Mg, xploring a weak (esN - 0.03) 
and a strong (esN - 1 feedback). The sizes of the systems are the 
same and their central gas density are respectively Pc.gas = 0.025 
and Pcgas - 0.066 me /cm-'. 

The initial number of particles went from 2'^ = 8'092 to 
222 _ 4' 2 94' 304, corresponding to a mass resolution of about 
10^ Mq and 20 M©, respectively. This is equivalent to changing 
the mass resolution by a factor 512, and the spatial resolution by 
a factor 8. We start by discussing the case of the weak feedback. 
The main properties of the simulations are summarized in Tab.[T} 
for all, Ci, = 0.05 and /vii - 0.25. Fig.|4]displays the evolution of 
the star formation rates and the corresponding increases in stel- 
lar mass. The amount of stars formed results from the complex 
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Table 1. Properties of the runs performed for the convergence 
tests. Mtot is the initial total mass of the system, Pc.gas is the ini- 
tial central gas density, is the initial total number of particles 
(gas and dark matter) and eg is the gravitational softening. These 
models were run with = 0.05 and esN = 0.03 





Pcgas 


N 






eg 


[lO'* Mo] 


[mn/cm^] 






[Mo] 


[pc] 


9.5 


0.071 




= 8' 192 


3.16 X lO-* 


100 






2l6 


= 65'536 


4.11 X 10' 


50 






2'" 


= 524'288 


5.00 X 10^ 


25 






222 


= 4'194'304 


6.17 X 10' 


12.5 


3.0 


0.022 


2U 


= 8' 192 


1 X 10* 


100 








= 65'536 


1.3 X 10' 


50 






2i9 


= 524'288 


1.58 X 10^ 


25 






2^2 


= 4'194'304 


1.95 X 10' 


12.5 



interplay between the gas cooling and feedback heating and thus 
traces well all the implemented physical processes. Fig.|4]shows 
that changing the resolution keep the peaks of star formation at 
the same epochs. High resolutions models form two to five times 
more stars over a longer period than the low resolution ones. This 
is the consequence of an improved sampling of the high central 
density peaks, which translates into shorter cooling times. 

The impact of a change in numerical resolution is slightly 
larger for small systems than for massive ones. The two simu- 
lations with 4'194'304 particles are still running at the time of 
writing this paper hence their corresponding red curves are still 
incomplete. However, during the first 4 Gyr, the star formation 
history of the 9.5 x lO'^IVIo model is similar for the - 2'^ and 
N = 2^^ resolutions. In the case of the smaller mass system, the 
difference in stellar mass between the = 2'^ and N - 2^^ mod- 
els is of the order of 15% during the 8 first Gyr. This increase in 
mass is nevertheless moderate and does not translate into differ- 
ent star formation histories, meaning different age or metallicity 
distributions and chemical abundance patterns. 

Fig. |5] show the stellar profiles, which depend directly on 
the dynamics. They show that a low spatial resolution leads to 
flattening of the profiles compared to higher resolution. 

Indeed, Fig. |6] illustrates how the metallicity distributions 
and abundance ratios are sensitive to the resolution. A low a res- 
olution (A^ < 2'^) tends to shift artificially the metallicity dis- 
tribution to lower values and only deliver an incomplete view of 
the chemical evolution. Satisfying convergence is again obtained 
for A^ = 2'^ particles, leading to [Fe/H] differences of less than 
0.1 dex. 

We performed the same tests for a strong feedback (esN = 1) 
and could verify that convergence is also satisfied, despite the 
fact that for these extreme feedback cases, inhibiting star forma- 
tion, the resulting low number of stars formed may lead to large 
noise as star formation is implemented as a stochastic process. 

In summary, our tests demonstrate the reliability of the nu- 
merical implementation of the physical processes described in 
Section |2] which are independent of the number of particles 
used. Above a resolution of A^ = 2'^ particles, the physical quan- 
tities are described with sufficient accuracy. The convergence 
arises from the balance between feedback and cooling; our sys- 
tems are self-regulated. In the rest of this paper, we use A^ - 2'^, 
a good trade-off between resolution and CPU time. 




s s 
Time [Gyr] 



Fig. 4. Star formation rate and stellar mass as a function of time 
for the two models of Table [T] Each line corresponds to a differ- 
ent resolution. 
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Fig. 5. Stellar density profiles for the models of mass 3 (left) and 
9.5 X IO^Mq (right) described in Table [l] with different resolu- 
tions, at f = 14 Gyr. 



4.3. Random number seed 

Whilst convergence of our results is obtained, nevertheless our 
systems could still be nonlinear at fixed resolution, due to the 
complexity of the physics. For example, changing the cooling 
impacts the star formation which in turn impacts the feedback 
and the cooling again, etc. The integration time of 14 Gyr is long 
with regards to the dynamical times of the galaxies, of the or- 
der of lOOMyr. Small perturbations might produce a large dis- 
persion in the final properties of the galaxies. In the following, 
we investigate this issue and evaluate the error bars attached to 
our outputs. Random numbers are used by the star formation al- 
gorithm and can act as perturbations. Therefore, we performed 
three different sets of simulations, in which only the random 
number seed was modified. The left part of Table |2] summarizes 
the parameters of our simulations. 
Fig. 



B.l 



displays the case of model #1111, which presents 
the largest deviations. They are anyway very small. Even after 
several dynamical times of evolution, the system is not strongly 
affected by the perturbations generated by the use of a different 
random number sequence. 
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Table 2. Properties of the models for which the random number seed was modified. 
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Fig. 6. [Fe/H] distributions and [Mg/Fe] vs [Fe/H] diagrams at 
/ - 14Gyr for the two models described in Table [T] The galaxy 
with a mass of 3 x 10** M© is shown on the left column, while the 
model with 9.5 x IO^Mq is displayed on the right column. 



Table |2] quantifies the relative variations in stellar mass, V- 
band luminosity, and final mode metallicity. These variations are 
estimated as the ratio between the absolute difference between 
the minimum and the maximum of the measured quantity di- 
vided by its mean value. The final stellar masses deviate by ~ 6% 
at most and the metallicities by 11%, while the luminosities can 
vary by 26%, due to the strong dependence of the luminosity on 
the stellar age. These values can be interpreted as the error bars 
intrinsic to the models. They must be kept in mind when drawing 
comparison between models and observations. 



5. Parametrization 

Physical processes involved in galaxy formation, such as star for- 
mation, feedback or cooling occur typically on lengths and time 
scales much smaller than the ones resolved by current numerical 
simulations. As described in Section|2] such sub-grid physics is 
treated through phenomenological prescriptions and requires the 
introduction of a set of parameters. Table [3] compiles the list of 
our model parameters. In the following, we explore their role. 
Unless specified otherwise, the discussions are based on the full 
set of 393 simulations. 



5.1. Feedback efficiency 

The feedback efficiency, esN, is the fraction of energy released 
by the explosion of supernovae, which is effectively deposited 
into the interstellar medium. The remaining fraction is assumed 
to be radiated away, without impacting the system. Our feed- 
back prescription is thermal (see Section 2.4 1, hence it modifies 
the temperature, the pressure and consequently the density of 
the gas surrounding the exploding stellar particles and in turn 
impacts the eligibility of the gas particles to further form stars. 
Fig. 7] displays the final mode stellar metallicity [Fe/H] and V- 
band luminosity (L,,) as a function of the initial total mass of the 
system. 

If 100% of the supernova energy is injected in the ISIVI 
(esN = 1, open circles in Fig|7]i, the final mode metallicity of 
the system is always below the values measured for the Local 
Group dSphs and UFDs (red vertical line). Similarly the lumi- 
nosities are generally very low. For such high values of esN, the 
feedback energy is maximal and strongly heats and dilutes the 
ISM pushing the gas particles away from the criteria for star for- 
mation, which is stopped. The subsequent increase of the cool- 
ing time, mainly driven by the decrease of the gas density, is so 
strong that it inhibits any star formation during the rest of the 
galaxy evolution. 

Clearly esN must be decreased below 1 in order to reach 
metallicities and luminosities compatible with the observations. 
Fig.jTjqualifies esN in the range 0.1 to 0.01. Besides, it also shows 
that the initial total mass of the systems drives the final luminos- 
ity and metallicities. esN primarily influences the least massive 
systems. At very low esN the cooling of the gas is no longer 
counterbalanced by the supernovae feedback: stars can form effi- 
ciently and longer, leading to high metallicities. Fig. B.2 and B.3 
illustrate also how the ISM becomes more homogeneous with 
decreasing esN. diminishing the dispersion in stellar abundance 
ratios. To further constrain esN, keeping its value fixed for all 
galaxies, requires to consider the full set of dSph properties. Our 
generic models (see Section 7.2 1 are obtained with esN = 0.03. 



5.2. Threshoid density for star formation psfr 

The threshold density for star formation, psfr, has been intro- 
duced in a cosmological context by Summers ( 1984| l in order 
to avoid the formation of stars in low density regions. The value 
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Table 3. List of model parameters. The reference values indicate our final set of adopted parameters based on the four generic 
models of Local Group dSphs. The last column indicates whether or not the parameters were varied (see Section|5]and Section|6]l, 
in which case we give their range of values. 



context 



quantity 


symbol 


baryonic fraction 


/b 


dark halo core radius 


rc 


virial fraction 


/vir 


total mass 




gas central density 


Pc tot 


initial outer radius 


'"max 


softening 


fg 


number of neighbors 




coolmg lunction 


A 


parameter 


c* 


critical density 


Psfr 


maximal number of stellar 




particles per gas particle 


A'* 


gas temperature 




initial mass function 


IMF 


yields 




supernova energy 


-EsN 


efficiency 


fSN 


adiabatic period 





reference values 



variation 



initial conditions 



gravity 
sph 

cooling 

star formation 



stellar evolution 



thermal feedback 



0.15 
1 kpc 
0.5 



50 pc 
50 

Sutherland & Dopita|(T993| ; 



Maio et al. 



(20071 



0.05 

0.1 mn/cm^ 

4 

3 X lO'' K 
Kroupa (■2001j_ 



Iwamoto et al.|^l'j99 l; Isujimoto et aL| ( |1995f 



10^' ergs 
0.03 
5Myr 



no 
no 
no 

[1 - 9x 10** Mo] 
[0.005 - 0.063 mn/cm^ ] 
[1.91 - 24.4 kpc] 

[adaptive] 

[50, 100] 



[0.01 - 1] 
[0.05- lOOmH/cm^] 

[1-15] 
no 

[Kroupa,Salpeter] 
no 

no 
[0.01 -1] 
[2 - 30 Myr] 
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Fig. 7. Final [Fe/H] and luminosity Ly as a function of the total 
initial mass of the system. Our full set of simulations is shown. 
The colors code the supernova feedback efficiency, between 0.01 
and L In each panel, a red line indicate the range of metallicity 
or luminosity covered by the Local Group dSphs. 



duced by |Katz et al.| ( |1996l l and corresponds to the mean density 
of the warm neutral medium of our galaxy. Since psfi fixes the 
local density at which stars are allowed to form, it may as well 
influence the global galaxy star formation history. 

The importance of this threshold density has been discussed 
in the context of disk galaxies ( [Tasker & Bryan|[2006[ |2008| l. 
Saitoh et al. (20081) concluded that only models using a high 
threshold density (psiv - lOOmH/cm"*) are able to reproduce the 
complex, inhomogeneous and multiphase structure of the ISM. 
Hereafter, we investigate the consequences of varying the critical 
density in the context of dSph galaxies. 

For c* = 0.05 and esN = 0.05, we varied psfr from 0.05 to 
100 mn/cm-' for two different initial masses, 3.5 and 7x10** Mq, 
within r^jax = 8 kpc and the initial central gas densities of 
0.025 and 0.053 me/cm^, respectively. A zoom on the first few 
Gyr evolution of the star formation rate and stellar mass of the 
3.5 X 10^ Mq model is provided in Fig. [s] It shows that the in- 
crease of Psfr imposes an increasing delay for the onset of the 
star formation, because the gas needs to be gradually denser. 
Moreover, one passes from a regime of large and moderately 
intense peaks of star formation to a series of higher intensity and 
higher frequency episodes, creating a larger number of stars, al- 
though varying psg by a factor 1000 does not change the final 
stellar mass by more than a factor 2.5. Drastic changes in the 
metallicity distributions and abundance ratios can be witnessed 



in Fig. B.8 



Alimi 



of 0.1 mn/cm^, that is widely used in the literature (e.g.,^ 

[eraL]|20U3l [Stinson et al.l2006l [Valcke et al.|[2008l l, was intro 



Following the course of the formation of the galaxies, it 
stands out that for low ps^ (< I mu/crn^), the first gener- 
ations of stars form within a relatively large region (radius 
~ 300 pc), roughly corresponding to the SPH sphere (radius 
~ 250 pc) in these low density medium. Metals and energy feed- 
back are homogeneously distributed. For larger star formation 
density thresholds, the evolution of the systems is quite different. 
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Time [Gyr] 



Fig. 8. The effect of varying psfr. This figure illustrates the case 
of a 3.5 X IO^Mq mass system, with an initial gas density of 
0.025 mn/cm^. The upper panel shows the evolution of the star 
formation rate with time, while the lower panel displays how 
mass increases. 



Because the density of the gas is higher, the size of the star form- 
ing region is smaller, with typical radii of the order of ~ 30 pc 
for lOOmn/cm^, again corresponding roughly to the SPH sphere 
radius (~ 15 pc). The mechanism of metals and energy feedback 
being linked to the SPH formalism, metals are distributed very 
locally. The release of the energy of the supernovae explosions 
creates bubbles. Those inflate and push back the surrounding 
gas, which significantly increases its density. The rims of the 
gas bubbles are dense enough to cool quickly. They become un- 
stable and form stars. In this way, star formation induces star 
formation. This explains the increase of stellar mass for large 
Psfr, observed in the bottom of Fig. |8] 

The metallicity distributions and abundances ratios gener- 
ated by Psfr > 1 me /cm-' look incompatible with the observed 
galaxy properties and are consequently discarded. We thus keep 
the value of 0.1 which best fit the abundances constraints. 



5.3. Star formation parameter c-i. 

The star formation parameter c,^ is a dimension-less parameter 
that controls the amount of stars formed from the gas. In the con- 
text of spiral galaxies, c* may be calibrated using a Kennicutt- 
Schmidt law ( Schmidt|1959[|Kennicutt 1998 ). Unfortunately, as 
dSphs are devoid of gas, no similar relation exists for them. The 
impact of the variation of has been discussed by |Stinson et al.| 
( |2006| l in the context of their feedback blast wave model. Their 
mean star formation rate varies by only a factor ~ 1 .2 when c* 
goes from 0.05 to 1 . 

We performed a similar study for our pressure supported sys- 
tems, varying c* by two dex, from 0.01 to 1. We considered two 
models with 3.5 x 10**Mo (pcgas = 0.025 mH/cm\ esN = 0.05) 
and 7 x IO'^Mq (pe,gas = 0.053 mH/cm^ esN = 0.05). The com- 
parison of the galaxy final state is provided in Fig.|B.6|and|B.7| 
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Fig. 9. Star formation rate and stellar mass as a function of time 
for models of mass 3.5 x lO^M© (pcgas = 0.025 mn/cm-') with 
different c*. 



Fig. |9] offers a zoom into the first few Gyr of evolution for the 
3.5 X 10*^ M0 model. Independent of the mass of the dSph, large 
results in a very discontinuous star formation history. This is 
a direct consequence of the coupling between star formation and 
feedback. In the extreme case of c* = 1 and Mtot = 3.5 x IO'^Mq, 
99% of the stars are formed during a first burst. This dramati- 
cally impacts the metallicity distribution and metal abundances 
(see Fig. 



B.61 



Metal-poor stars ([Fe/H] < -2.5) with high 
[Mg/Fe] ~ 1 ratio are dominant and make these systems unreal- 
istic. The long period of quiescence following the first and domi- 
nant peak results from an important physical process, which will 
be discussed in Section fTTI 

On the contrary, at lower c* the galaxy experiences a more 
continuous star formation history, resulting from a better self- 
regulation between cooling and heating. The moderate intensity 
of the initial bursts allows the gas to cool. This gives rise to the 
presence of intermediate age stellar populations and higher mean 
metallicity. 

The full evolution of the stellar mass with time is displayed 
in Fig. B.6 and B.7 Differences in the final stellar masses remain 
smaller than 50%. Increasing slightly decreases the quantity 
of stars formed. This counter intuitive result is also found by 
Stinson et al. (2006 1 at large c*. Indeed, a model with a higher 
forms quickly a large quantity of stars. The feedback energy 
due to the SNell explosions is released nearly instantaneously. 
The gas is blown out and the star formation is quenched. At later 
times, the large number of SNela continue to prevent star forma- 
tion, as discussed in more details in Section [7T| 

5.4. Number of stellar particles created from a gas 
particle 

The mass of the stellar particles m^, is linked to the mass of the 
gas particles nig via the parameter A^* = mg/m*. It can also be 
seen as the maximum number of stellar particles that can be ere- 
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ated per gas particle in absence of accretion of gas. In princi- 
ple, Eq. |4] regulates the impact of the resolution in mass of the 
simulations, leading to the same integrated star formation rate. 
Nevertheless, this is not fully the case once cooling and feedback 
are taken into account. Augmenting A^* increases the sampling 
of the star formation. A small number of big particles are re- 
placed by a larger number of smaller ones, distributed along the 
same period of time, and for an equal final stellar mass. The price 
to pay is obviously a large increase of the CPU time needed for 
a simulation. According tolSpringel & Hemquist ( 2003| l, a good 
compromise between time resolution and CPU saving is Ni, = 4. 

We varied A^* between 1 and 15 for a 5 x IO^Mq (rmax - 
8kpc, Pcgas = 0.037 mu/c nv', = 0.025, esN = 0.02). The re- 



sults are displayed in Fig. B.IO Two effects are observed. The 



final mode stellar metallicity of the systems decreases within 
larger A^*, by a factor two passing from A^*=l to 15. A large 
number of small stellar particles regularly injects feedback en- 
ergy in the interstellar medium, preventing it to cool and be el- 
igible to star formation again. Conversely, as the probability of 
forming stars for large A^* is small, the system is slightly cooler 
and denser when the first star form, leading to a higher star for- 
mation peak. Whilst the mean metallicity of the systems are 
similar within 0.3 dex, a clear trend is observed in the abun- 
dance ratios. The dispersion in [Mg/Fe] is significantly reduced 
for higher A^* (high chemical mixing), again owing to the fre- 
quent feedback of the numerous stellar particles, starting early 
in the evolution of the galaxies. 

Small abundances dispersions at low metallicities are the 
rule in observed galaxies. This favors high A^,^. The dispersion 
does not increase substantially between Ni, - 4 and A^* = 15 
(see Fig. |B.10] l, therefore saving CPU time, we set it to 4. 

5.5. Supernova adiabatic period t^d 

The adiabatic period f^d fixes the period of time during which 
the cooling of the gas particles, that have received the thermal 
supernova feedback, is stopped. As discussed in Section 2.4 this 



period improves the feedback efficiency by avoiding the instan- 
taneous loss of energy by radiative cooling, fad must be much 
larger than the dynamical time steps, which is at most 0.2 Myr 
in our case as imposed by the resolution of the chemical evolu- 
tion (see Section [23] l, in order to ensure the convergence of the 
simulation. It is often set to 30 Myr, the shortest lifetime of the 
least massive SNII progenitor (8 Mq). Stinson et al. ( 2006^ pro- 
posed a new feedback recipe based on an analytical treatment of 
supernova blast-waves. In this model, fad depends on the local 
gas density and pressure. 

We explored the effect of fad, between 2 and 30 Myr, for two 
models, Mtot = 3.5 (r^ax = 8kpc, Pcgas = 0.025 me/cm-') and 
7 X 10** Mq (r^ax - 8 kpc, pc.gas - 0.051 mn/cm-*). Both models 
have c* = 0.05 and esN - 0.05. The results are displayed in 



Fig. B.4 and B.5 Increasing fad decreases the final stellar mass in 
both cases. This is the consequence of a lower cooling. However, 
the effect remains weak. After 14Gyr, the amount of the stellar 
mass formed is diminished by a factor of 1.3 to 1.6. This is not 
sufficient to change substantially the chemical properties of the 
final dwarf. In the following, we fix fad to 5 Myr. 

5.6. Number of particles in a softening radius A^ngb 



As described in Section [23] and \2A\ the stellar ejecta as well as 
the feedback energy are spread in the SPH sphere. Increasing the 
size of the SPH sphere, i.e., increasing the number of neighbors 
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Fig. 10. The difference in stellar mass growth with time 
when changing of IMF. Two examples are shown, at 2 x 
lO^Mo (pcgas = 0.015 mn/cm^) and 8 x IO'^Mq (pc.gas = 
0.053 mn/cm^). 

when computing the local physical quantities, allow to distribute 
metals at larger radii. 

In galaxy cluster simulations, [Tornatore et alT] (|2007) varied 
A^ngb from 16 to 128. They reported only marginal differences, 
although increasing the number of neighbors produced some- 
what higher star formation rates, presumably due to the largest 
number of particles affected by metal-line cooling. 

We performed a similar test with a 4 x 10^ M© (r^ax - 8 kpc, 
Pccas - 0.029 mn/cm-*) simulation, and changing A^ngb from 50 
to'loO 



The final state of the systems are display in Fig. B.13 
The difference between the two models is very small, in agree 
ment with |Tornatore et aL] ( |2007 )l. However, we reach an oppo 
site conclusion. When doubling A^ngb> the final stellar content is 
slightly decreased (factor 1.33). This induces a slight decrease 
of the mean metallicity (< 0.1 dex). The reason is probably 
linked to the different treatment of the feedback in the two stud- 



ies. 



Tornatore et al.^2007 )'s model is based on the wind scheme 



of [Springel & Hemquist ( 2003 | l. 



5.7. The IMF 

We finally checked the impact of the choice of IMF on our 
results. For this, we used the classical Salpeter IMF (Salpeter 
1955 1, defined by a constant slope of -2.35 . The 'Kroupa 
( |2001| l's IMF contains less stars with masses smaller that 0.2 Mq 
than the Salpeter's one. As a consequence of normalization, this 
means more stars in the mass range of supernovae (> 8 M©) by 
a factor 1.5, and thus more feedback energy. Fig. [T0| illustrates 
that less stars are formed with a Kroupa IMF than with a Salpeter 
one. 

This weak difference in the final stellar mass does not impact 
substantially neither the metallicity distributions nor the [Mg/Fe] 



abundance ratios (see Fig. B. 1 1 and B. 12 1 



6. Mass and central density 

In the previous Section, we explored the role of the parameters 
inherent to GEAR. We now discuss the effect of the initial total 
mass of the galaxy and of the initial central gas density, which 
are linked to the initial conditions. 
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|Revaz et~aL ( |2009| l found the galaxy total initial mass esti- 
mated within a fixed radius to be a dominant parameter, increas- 
ing both the final metallicity and the luminosity of the galaxies. 
We confirm this result as illustrated in Fig.|7] which can be com- 
pared with the Fig. 5 of Revaz et al. (2009 ) . 

As described in Section 1 3. 1| the dark matter profiles issued 
from the ACDM simulations exhibit a small dispersion in their 
central regions. This deserves investigation. Indeed, the initial 
gas profile follows the dark matter one in our models, and the gas 
density determines the cooling time. Moreover, the star forma- 
tion criterion involves a threshold in gas density. All these facts 
motivated us to examine the role played by the central density 
in addition to that of the total mass. In a previous publication, 
ICarraro et al.| ( |2001| l had observed different star formation his- 
tories for objects of the same total mass, when their the collapse 
phase started at different initial densities. 

At a given central density (pc.gas), the total mass (Mtot) of 
the system can be varied by increasing/decreasing its size with 
Tmax- Similarly, for a given total mass, the central density can be 
varied increasing/decreasing its size with r^ax- Fig. 11 presents 
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the variation of the final mode metallicity as a function of the 
initial mass and the central gas density. We present 40 simu- 
lations - all with = 0.05 and csn - 0.03 -, where Pc.gas 
is varied from 0.007 up to 0.063 mn/cm^ and Mtot from 1 to 
9x10** Mq. Consequently, the outer galaxy radius varies between 
1.8 to 25kpc. 

For the range of metallicities covered by the Local Group 
dSphs, an increase of the mass by a factor 10, at fixed central 
density, increases the final metallicity of the system by only 
a few tenth of dex. The largest variation is obtained for very 
small systems: from 10^ Mq to 3 x 10** Mg, the mode metallic- 
ity can increase by 0.4 dex at fixed density; above 3 x 10^ Mq, 
variations soften until they are essentially indiscernible. Above 
3 X 10** Mq, the additional amount of matter provided by the ex- 
tension of the system (increase of r^ax) does not affect much its 
final chemical properties (metallicity distribution and [Mg/Fe] 
vs [Fe/H]). It acts primarily on enhancing the final total stellar 
mass, while it slightly modifies the age distribution of the stellar 
population. This latter consequence is mainly seen at low densi- 
ties ( pcgas ~ 0.015 mH/cm^*). These points are also illustrated in 
Figs. B.14 and B.15 with different and esN than in Fig. 11 



Conversely, a variation by a factor 10 in the central gas den- 
sity can vary the final metallicity by more than a dex, making 
it the primary driving parameter. The cooling time is shorter 
for larger densities, enhancing the star formation, and there- 
fore resulting in more luminous and metal-rich systems. At high 
enough central densities, all models, even the least massive ones, 
may experience a strong initial burst, as seen in Fig. B.16I In the 
0.015 to 0.037 mu/crn^, the star formation 



range p, 



c.gas 



history 

can be fully modified. Example is provided in Fig. |B.17] where 
the central density is varied by a factor 4, passing from a low and 
continuous star formation to an initial strong burst. 



7. The different observed regimes of star formation 



Fig 12 presents in [Fe/H] and M/L vs luminosity plots the full 
series of simulations, which are globally compatible with the ob- 
servations. This represents a total number of 160 simulations 
with in the range 0.05 - 0.1, in 0.03 - 0.05, Pc.gas in 
0.005 - 0.06mH/cm3 and Mtot in 1 - 9 x 10'^ Mq. While the 
reproduction of these global relations is a necessary step, it is 
however not sufficient to ensure that the star formation history 
of the galaxies are correct. 
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Fig. 11. Final mode [Fe/H] as a function of the total initial mass 
of the system. The size of the circles scales with the initial size 
(radius) of the system, as featured at the bottom right of the 
figure. The color code the initial gas density. For all models 
= 0.05 and esN=0.03. 



Among these simulations, we selected a set of 4 models 
best reproducing the observations. We limited as much as possi- 
ble variations in esN and c*. Our selection was based on com- 
parison with the observed galaxy metallicity distribution, the 
V-band luminosities, the abundance ratios, and the stellar age 
distributions. Observed luminosities were taken from Wa lker! 
\20 09b|l. T he metallicity distributions were retrieved from 



etaT 



Batta gUa et al.| ( |2011.) for Sextans , from Battaglia et al.| ( [2006l l 



Tols toy et al.| p004| ), and Hel mi et al.| ( |2006| l, fOT Fornax, 



Sculptor, Carina, respectively, with the new CaT calibration of 
Starkenburg et al. (2010 1. The abundance ratios were taken from 
iToIstoy et al.|(|2 009'),'L etarte et al.| ( [20T0l i, |Shetrone"etaL] ( p003) , 
and Koch et al. (2008 ). The stellar age distributions were taken 
from |Smecker-Hane eiH] (|1996|l; [Hurley- Keller et al. (1998|l 
(Carina), Babusiaux et al. (|2005 ); Shetrone et al.|()2003| ); (Tolstoy| 



| eraL](p003]i (Sculpto r), |CoIeman & de Jong| ( |2608"] r(Fomax) , 
and |Lee et al. ( 2003 | l (Sextans). A summary of the mean ob- 
served properties of the four local group dSphs is given in Tab. [5] 



It also includes the dSph tidal radii from Irwin & Hatzidimitriou 
( jl995j) as well as the velocity dispersions from Walker et al. 
P009a^ . 



Fig|T3] displays our final selection. These models are shown 
with green squares in Fig. [T2j they fall close to the observed 
M/L, [Fe/H] and Ly, within a factor 2 for each quantity. Table 
|5]provides the model parameters and outputs, is the radius en- 
compassing 90% of the V-band luminosity. Mgas, Mstar, Mhaio and 
M/L are computed inside rt. [Fe/H] is calculated as the mode of 
the galaxy metallicity. The stellar and dark matter velocity dis- 
persions, (Ti, and ctdm, are measured along the line of sight, in- 
side the effective radii (half light radius). We differentiated the 
dark matter and the stellar velocity dispersions, because the for- 
mer best traces the total mass of the system, while the observa- 
tions gather the latter. 
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7.1. Distinct fates for tfie SNell and SNela feedbacks 

In all systems, the star formation history results from the bal- 
ance between cooling and heating. Both SNell and SNela play 
important roles, despite their relative different numbers and con- 
sequently different integrated energy feedback (a factor 10 be- 
tween the two). 

Looking closely at the course of the galaxy evolution, one 
sees that star formation is ignited in the galaxy central parts. 
Massive stars explode rapidly as SNell and their injections of 
energy counteract the increase of the central gas density, leading 
to a flat gas profile within ~ 1 kpc. The strong and quick energy 
injection of SNell explains why the rapid initial rise in intensity 
of the star formation rates (see first panels in Fig. 13 i is soon 



followed by a sharp decrease: the feedback energy acts against 
star formation. 

Meanwhile, lower mass stars are progressively redistributed 
along a steep profile. Indeed, their velocity dispersion at birth is 
small, of the order of 5 km/s, as inherited from the turbulent gas 
motion. Therefore their ensemble contracts because, being de- 
coupled from the gas, it is no longer supported by pressure. The 
number SNela increases with time until they become the domi- 
nant source of heating in the central regions (0.2 pc) preventing 
the gas to cool and condense again, despite a very short cooling 
time. The dominance of the heating by SNela over the SNell's 
one is only possible because the former are more concentrated 
than the latter This causes the few Gyr-long periods of very low 
level (or absence) of star formation after the initial 2-4 Gyr. 
Once the bulk of SNela explosions has passed, star formation is 
possible again at higher rates. 

These phenomena have been explicitly witnessed in a 
Sculptor-like simulation in which the energy released by SNela 
was artificially set to zero. Star formation was then maintained 
to a substantial level of 0.005 Mg/yr, corresponding to half of 
the amplitude of the first burst, where star formation had been 
quenched in the generic model. 

7.2. Generic models 

Overall, our models mean properties are in excellent agreement 
with the observations. This suggests that a sequence of dSphs 
can be reproduced by varying very few parameters. Carina, 
Sextans, Sculptor and Fornax form an increasing sequence of 
total masses and initial central densities. All systems display a 
first ~ 2 Gyr burst. Depending on the mass, the star formation 
history can be then considerably reduced, or even quenched, in 
the case of small systems. esN is fixed to 0.03; c^, - 0.05 except 
for Carina. Clearly, as discussed in Section 5.3 increasing 
favors discontinuous star formation histories. This small differ- 
ence may be the sign of some external factor in the formation of 
the galaxy, which was not taken into account in our models. 

The sequence in mass and central gas density can be under- 
stood by the following considerations, at fixed c^, {~ 0.05) and 
esN (~ 0.03): 

- At high mass, i.e., above ~ 3 x 10^ Mq. The evolution of 
the systems is quite uniform. The cooling always dominates 
the feedback and the virial temperature reaches lO'* K around 
which the radiative cooling is very strong. Hence, the star 
formation is never completely quenched. 

- The diversity of evolutionary paths is larger for lower mass 
systems. Those who initially form a large quantity of stars 
compared to their total mass, do release a large SNela-driven 
feedback energy, which is able to counterbalance their cool- 
ing. The level at which star formation is quenched depends 



dSph 


Lv 


<[Fe/H]> 


MIL 




cr 




[10* Le] 






[kpc] 


[km/s] 


Fornax 


14 


-1.17 


12 


2.08 


11.7 


Sculptor 


1.4 


-1.96 


158 


1.33 


9.2 


Sextans 


0.41 


-2.26 


19 


3.10 


7.9 


Carina 


0.24 


-1.86 


88 


0.58 


6.6 



Table 5. Some of the observed properties of the Local Group 
dSphs Fornax, Sculptor, Sextans, and Carina, which are used 
in the present study . The luminosity comes from Walker 
al? (2009b). The mean metallicities were calculated from 



et 



BattagHa et al7|( 2011 
(Tolstoy et aL| 2004 



for Se xtans, from (Battaglia et al.|20061 l, 



and (Helmi et 



2006 1, for Fornax, 



Sculptor, Carina, respectively, with the new CaT calibration of 
Star kenburg et al. ( 2010| l. The M/L ratios we r e computed us- 
ing the masses derived by Walker et al. (2007 1; Battaglia et al. 
(2008|;|Kleyna et al.| (|2004| l. The tidal and core radii are taken 



from Irwin & Hatzidimitriou ( 1995 i and the velocity dispersions 
from .Walker et al., ( ,2009ar "^ 



on the ratio between the cooling and the energy released by 
SNela, following the process described in Section 7.1 This 



ratio reflects the level at which the galaxy has used its full 
capacity to form stars during the first Gyr, i.e, the fraction of 
gas transformed into stars, as determined by the galaxy mass 
and initial central density. Low mass with low central gas 
density systems have a rather long cooling time (> 1 Gyr). 
The SNell are sufficient to completely quench momentar- 
ily the star formation, which needs more than 1 Gyr to be 
reignited, and proceeds further episodically. 

The four models reproducing the observed properties of 
Carina, Sextans, Sculptor and Fornax follow these rules. An ad- 
ditional feature is that the choice esN (primarily) and c* was 
definitely driven by Fornax. Its luminosity and mass require a 
vigorous and continuous star formation, in turn forcing both pa- 
rameters to be low. 

Just as in Revaz et al. ( |2009[ l, our models retain large masses 
of gas, from 10" to 10' M© and in inverse proportion of the 
galaxy stellar mass, i.e., low stellar mass systems have the largest 
fraction of gas at 14 Gyr. Moreover, both Sextans and Sculptor 
models need to have their star formation artificially stopped, re- 
spectively at ~ 5 and ~ 9 Gyr. The low esN and c* imposed 
by Fornax are inefficient in quenching it either by the means of 
thermal SNe feedback or even galactic winds. This strongly sup- 
ports stripping of the gas by external gravitational and/or hydro- 

( |2n06| ). 



dynamical interactions (see for example, Mayer et al. 



even though the bulk of the galaxy formation could be driven by 
their initial conditions. 



8. Spatial distribution of the stellar populations 

The presence of metallicity gradients in dSphs is debated. 
Evidence for stellar population segregation is reported in 
Sculptor (Tolstoy et al."2004), Fornax (Battag fia et~aL] [2006) 1, 
and Sextans (Battaglia et al. 2011 ). In these galaxies, the most 



metal-rich stars are concentrated in the galaxy central regions 
while the metal-poor ones are essentially found at every radius. 
Conversely no radial change in the metallicity is found in Leol 



dKoch et al.|2007bl|Bosler et al.|2007 l. 'Gullieuszik et al. 2009|l, 
LeoII jKoch et al |l2007a| l, and CVnl (Ural et al.,2010 ). Kirby] 



et al. (2011 1 reported low metallicity gradients of at most -0.21 



dex per core radius in LeoII 
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Fornax 


1335 


7 


0.059 


7.1 


0.05 


0.03 




13.9 


-1.01 


1.98 


9.4 


15.0 


2.4 


1.35 


8.80 


Sculptor 


1324 


5 


0.029 


9.6 


0.05 


0.03 


9.1 


1.50 


-1.75 


2.93 


6.4 


11.7 


1.9 


0.34 


4.45 


Sextans 


1316 


3 


0.022 


8.0 


0.05 


0.03 


4.7 


0.37 


-2.09 


1.58 


4.2 


9.7 


0.5 


0.07 


1.04 


Carina 


1281 


1 


0.022 


3.5 


0.1 


0.03 




0.24 


-1.93 


0.76 


3.1 


7.2 


0.2 


0.02 


0.63 



Table 4. Properties of the four best generic models for the Local Group dSphs, Fornax, Sculptor, Sextans and Carina. Left-ward 
of the vertical line are the physical input parameters of the models. Right-ward of the same line are the final outputs of the models. 
They are calculated within the radius contairung 90% of the galaxy total light. 
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Fig. 13. Properties of the four selected models representing Fornax, Sculptor, Sextans and Carina. From left to right, the plots display, 
the star formation rate and the evolution of the stellar mass, the normalized stellar age distribution, together with the evolution of 
[Fe/H]. The final [Fe/H] distribution and ratio [Mg/Fe] as a function of [Fe/H]. The grey region in the stellar age distribution of 
Sculptor and Sextans correspond to the stellar population which would have been created if we would not have truncated the star 
formation. These populations are absent from the metallicity distributions and abundance ratios. The stellar V-band luminosities 
computed are 14 Gyr. 
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Fig. 12. Global relations. The galaxy mass-to-light ratios, M/Ly 
(top panel) and the mode of their final metallicity distributions 
(bottom panel), versus their V-band luminosities Ly Each circle 
stands for one of our 160 models, with c* and esN restricted to 
the intervals 0.05 - 0.1 and 0.03 - 0.05, respectively. The cen- 
tral gas density and total mass range is 0.005 - 0.06 me/cm-' 
and 1 - 9 X 10^ Mq. Colors code the initial central gas den- 
sity. The red large squares represent the Local Group classical 
dSphs, while the smallest red squares stand for the ultra faint 
dwarfs. The dashed lines correspond to our model best fits. The 
green squares indicate the position of our selected generic mod- 
els, from left to right. Carina, Sextans, Sculptor and Fornax. 



Fig. 14 presents the radial [Fe/H] density maps for our four 
generic models. They show no evidence for any radial change 
in metallicity distribution. This is a common feature to all our 
models. In the following, we examine the gas motion leading to 
the final structure of our model galaxies and consider two of our 
very high resolution simulations containing 4' 194' 304 particles 
with an initial total mass of 3 x 10*^ M© and 9.5 x 10^ Mq, from 
Table [U 

The 3 X 10^ Mq model is characterized by two strong 
episodes of star formation (see the top of Fig. |4]). The first one 
occurs between 0.4 and 1 .6 Gyr, the second one, slightly weaker, 
takes place between 2.1 Gyr and 3.2 Gyr. The physical processes 
described in the following apply to any other model with a small 
mass. Fig{T5] displays the radial [Fe/H] density maps of the gas 
during the first period of star formation. The metallicity gradient 
seen in the first snapshot at f = 0.8 Gyr is the consequence of 
the first generations of stars forming preferentially in the galaxy 
central high density regions, hence leading to higher chemical 
enrichment. This implies that for a short period of time, stellar 
evolution occurs on shorter time scales than the gas motion. The 
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Fig. 14. The radial distribution of the stellar [Fe/H] for our four 
selected generic models. Each panel is a 2D histogram of the 
number of particles weighted by their mass. The blue curve fol- 
lows the mean value of [Fe/H] computed at each radius bin. 



gradient persists at f = 1.4 Gyr, although being already shal- 
lower, but soon after it disappears. 
Fig. 



16 presents the physical mechanisms at play. As seen in 
the upper panel, the gas accumulate in the central regions during 
the period of intense star formation, not only metals but also the 
thermal energy released by the explosions SNell. Consequently, 
the gas is hotter and more tenuous than its surrounding (see snap- 
shot at f = 1 .4 Gyr) and forms a bubble. Very strong Archimedes 
forces act in the vicinity of the center of the galaxy potential 
well, making the position of the bubble unstable. Small motions 
induced by the local turbulence is sufficient to offset it and makes 
it quickly dragged outward. This phenomenon is seen twice: first 
between 1 .5 and 1 .7 Gyr, then between 1 .7 and 2 Gyr. The hot 
and metal-rich gas bubble (seen in red in Fig. 16 1 leaves the main 
body of the galaxy, reaching regions up to 10 kpc. Subsequent 
episodes of star formation would experience the same scenario, 
erasing any metallicity gradient in the gas in less than a Gyr. 
We note that the uplift of hot bubbles by Archimedes forces is a 
common mechanism in cooling flow clusters (Revaz et al. 2008| l. 

The evolution of the 9.5 x 10** Mq model is displayed at the 
bottom of Fig. 



16 



A gradient in metallicity is visible between 
1 . 1 and 1 .7 Gyr. The fundamental difference with the former ex- 
ample is that here stars efficiently form in a much wider region. 
Metals and feedback energy are then injected in a larger volume. 
Therefore, turbulent motions are sufficiently strong to erase the 
initial gradients. 

In summary, our high resolution simulations show that the 
hot gas motion have very short timescales of evolution, that are 
incompatible with the formation of stellar metallicity gradient. 
We note that our simulations do not include mixing in the in- 
terstellar medium. This would dilute even more any differential 
spatial metallicity distribution. Conversely, a multi-phase and 
multi-scale structure of the interstellar medium may prevent the 
formation and migration of hot bubbles. This probably deserves 



Y. Revaz and P. Jablonka: The Dynamical and Chemical Evolution of Dwarf Spheroidal Galaxies with GEAR 



15 



future investigation. As discussed in Sect. |5.1| the effect of our 
feedback is rather small. In principle, a stronger feedback would 
increase the turbulence of the ISM and also fight against metal- 
licity gradients. 



9. The stellar mass fractions in N-body simulation 
dark matter haloes 

The final baryonic fraction fbf is an interesting quantity, since it 
can be observationally derived in galaxies and is diff'erent from 
the initial cosmological one. Within the radius containing 90% 
of the stars, this fraction varies in our models from 25% to 35%. 
Considering the stellar baryonic mass, it ranges from 2% to 10%, 
similarly to the results of Valcke et al. ( 2008 1, who ran models 
of galaxies in isolation as well. 

Looking at Mstais/Mhaio, |Sawala et al.| ( [2011| ) compared the 
predictions of a series of N-body models of galaxies with the 
expected ratios calculated from the Sloan Digital Sky Survey 
(SDSS) for the stellar part and from the Milleniumll simulation 
for the pure dark matter content [Gno et al. 2010). The depen- 
dence of Mstars/A^haio on Mhalo was Originally performed for stel- 
lar masses between IO^^Mq and IO^'^Mq. It wa s then extrap- 



olated from 10'* '*Mo down to IO^'Mq. Doing so, [Sawala et al. 
( |201 1) found the fraction of stellar baryons in dark halo masses 
below IO'^Mq to be lower in the SDSS/Milleniumll predictions 
than in the galaxy models. Their interpretation was that galaxy 
models may produce too many stars. This conclusion deserves 
further consideration. 

Figure [17] moves the predictions of |Sawala et al.| ( |201 l| l 
into observable quantities. Guo et al. (2010rs relation between 



galaxy stellar and halo masses was turned into V-band luminosi- 
ties and velocity dispersions. Following I Sawala et al.| ( |201 l| l, 
we consider two different slopes for the faint-end of the stellar 



mass function, -1.15 (|Li & White|2009| and -1.58 (Baldry et al. 



[2008 ). Turning stellar masses into luminosities was done with a 
reasonable stellar Mstars/Lv ratio of 0.75 ( Flynn et al.|2006 l. As 
the observed V-band luminosities vary over several dex, uncer- 
tainties on Mstais/Lv impact hardly the results. The correspon- 
dence between the halo masses and their velocities was com- 
puted in two ways: first, we used the relation between the max- 
imal velocity and halo mass measured in the Aquarius project 
( Springel et al.|2008] ). Second, we assumed an NFW profile and 
solved numerically the Jeans equations in spherical coordinates 
following Eq. [TT| and considering extreme values, 5 and 30, for 
the profile concentration ( jMaccio et al.|2007] l. The velocity dis- 
persions extracted from the inner 1 kpc are computed by prop- 
erly taking into account projection effects. 

The comparison with the observations - dwarf galaxies 



(Mateo 1998 1, together with irregular, elliptical and spiral galax- 



ies from (Garrido et al. 
[Zeeuw et al.p002 



2002) 1 and the SAURON sample (jdej 



and a range of dififerent N-body simula- 



tions ( ISawala et al.|201 l|[Stinson et al.|2007|[20()9l[Valcke et al.| 



2008] [Govemato et al.||2010[ [Pelupessy et al.||2004| i call for a 
few comments: Noteworthily, the agreement between all mod- 
els of small mass systems on the one hand and their consistency 
with the observations on the other hand is quite remarkable, de- 
spite considerable variations in the assumptions, including the 
initial conditions, from cosmological to isolation. Where stellar 
masses were actually measured in the SDSS, the consistency of 
|Guo et alT] ( |20lO| )'s relation with the observations is very good in- 
deed. Conversely, at lower masses, the discrepancy with Sawala 
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Fig. 17. The relation between the galaxy velocity dispersions or 
maximal velocities and their V-band luminosities. Filled circles 
correspond to the results of this work and those of N-body sim- 
ulations. The observations of low-mass galaxies are displayed 
with open symbols. The gray area indicates how the Guo et al. 
( |2010| l's relation between stellar and dark halo masses, further 
extrapolated by |Sawala et al. |(|201 l|l, translates into luminosities 
and velocity dispersions for a minimum and maximum value 
of the NFW concentration parameter. Below a stellar mass of 
IO^-^Mq, two faint end slopes of stellar mass functions are con- 
sidered (see text), splitting in two diff'erent branches the grey 
area. The dashed black lines are obtained when considering the 
Aquarius project's relation between the maximal velocity and 
halo mass ( Springel et al.l[2008[l. The min imum stellar mass 



(lO'^ ^'Mo ) for the original |Guo et al. (2010 1's fit and its corre 



spondence in terms of range of velocities is indicated with black 
arrows. 



the dark matter halo ones by a factor 2 (see Table |4|, this is in- 
sufficient to fill in the gap between the SDSS/Milleniumll rela- 
tion and the observations. Similarly, varying the stellar Mstars/ Ly 
even by a factor 10 does help either. 

The velocity dispersion of the stellar component depends 
also on its density profile, which could well be different from 
the dark matter's one. Therefore we performed some tests as- 
suming different profiles for the stars and the dark haloes. For 
the stellar density profile, we took a truncated plummer sphere 
(e.g., |Battaglia et al.|2006| [2008 ), while the dark matter density 
profile was kept to the NFW form. We calculated the inner Ikpc 
velocity dispersions of both components; they agree to within a 
few percents. This is much lower than the mean factors observed 
in Fig. 17 



et al. (2011 1's extrapolation is very large. Whilst the observed 



velocities are derived from baryonic matter and could undermine 



In summary, the relation between Mjtars derived from ob- 
served luminosities, color, spectroscopy etc., and Mhaio from 
pure dark matter simulations should likely be revisited. The as- 
sumption under which the number density of galaxies and dark 
matter haloes match each other ( |Guo et al.|2010| l lead to incon- 
sistency with the observations at low stellar masses. 
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Fig. 15. Evolution of the stellar metallicity gradient with time for a 3 x 10** Mq model containing 4' 194'304 particles. As in Fig. 
each plot corresponds to a 2D histogram of the number of particles weighted by their mass. 
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10. Conclusions 

We described our new parallel Nbody/Tree-SPH code, GEAR. 
Starting from the public version of Gadget-2, we included the 
physics of baryons, i.e., metal-dependent gas cooling, star for- 
mation recipes, thermal feedback due to Type la and Type II 
supernovae explosions, and chemical evolution. GEAR has been 
conceived to sustain high spatial resolution together with de- 
tailed chemical diagnostics, and to follow the galaxy evolution 
over a full Hubble time, either in isolation or in a cosmological 
context. We qualified his performances with the case of dSph 
galaxies. 

> GEAR conserves the total energy of the systems to bet- 
ter than 5% over 14Gyr. Moreover, the code proved excellent 
convergence of the results with numerical resolution. We could 
show that for small galactic systems, such as dSphs, 2"' parti- 
cles offer a good trade-off between CPU time and resolution, 
catching the essence of the physics at play and allowing secure 
predictions. 

> Varying the initial random number seed, we estimated that 
our intrinsic errors were below 30% on the model stellar masses, 
the V-band luminosities and mode metallicities. 

> Hundreds of simulations were performed in order to under- 
stand and quantify the effect of the free parameters, such as the 
star formation parameter (c*), the star formation density thresh- 
old (psfr), the supernova efficiency (csn), the number of stellar 
particles formed from each gas particle (A^*), the adiabatic pe- 
riod (fad) attached to the thermal feedback,and the number of par- 
ticles in a softening radius (A^ngb)- The most sensitive parameter 
is definitely esN^ with an acceptable range of values between 0.03 
and 0.05 to reproduce the sequence of dSphs. These low values 
further imply that strong winds are incompatible with their ob- 
served metallicities. 

> We ran a 5 12-' ACDM cosmological simulation of pure 
dark matter in order to study the profiles of haloes with masses 
between 10** and lO'^M©, that are typical of dSphs. The expan- 
sion of the Universe leads to the formation of stable systems after 
z = 6. In turn, this justifies models of dSphs in a static Euclidean 
space, where the expansion of the universe is neglected. The 
physics of baryons that depends on the density in physical co- 
ordinates is correct. We showed that haloes experiencing only 
minor mergers since z = 6 do exist. Whilst at fixed mass, the den- 
sities of these haloes are very similar, still they exhibit a small 
dispersion (factor 3 to 4). 

> We confirmed that the total initial mass plays a primordial 
role in the evolution of the galaxies. However, we demonstrated 
that the initial central gas density is as crucial. Changing the 
mass by a factor 10 translates into a few tenth of a dex in the 



final mean metallicity, while a factor 10 in central density can 
increase [Fe/H] by more than 1 dex. We showed that the Local 
Group classical dSph could have their properties reproduced as 
a sequence of mass and density. 

> Differences in mass and initial gas central densities seen in 
the ACDM simulation lead to a variety of star formation histo- 
ries explaining the diversity of the chemical properties observed 
in dSphs. In massive and dense systems, the cooling dominates 
the feedback, and stars are form continuously, leading to lumi- 
nous and metal-rich galaxies. At lower masses, the variety of 
star formation results from the subtle balance between the cool- 
ing and supernovae feedback of both SNela and SNell. Gas is 
still present in our model galaxies with masses of about 10* to 
10' Mq after 14Gyr. This strongly supports the existence exter- 
nal processes like tidal or ram pressure stripping, which we did 
not include in this study. 

> Because SNela and SNell have different spatial distribu- 
tions, they act in distinct ways in the heating of the systems. In 
particular, we showed that the few Gyr-long quiescent periods of 
star formation are due to the explosion of the centrally concen- 
trated SNela. 

> We investigated the relationship between the stellar mass 
of galaxies and their parent dark matter haloes. The agreement 
between all models of small mass systems on the one hand and 
with the observations on the other hand is quite remarkable, de- 
spite considerable variations in the assumptions, including the 
initial conditions, from cosmological to isolation. They are con- 
sistent with the observations. 

> Despite the fact that the new stars are preferentially formed 
in the galaxy central regions, none of our models displays any 
segregation in stellar population, such as the metal-rich being 
more centrally concentrated than the less metallic ones, as re- 
ported in some observations. We showed how turbulence in mas- 
sive systems and uplifted hot metal-rich bubbles in less massive 
ones do erase the initial metallicity gradients seen in the gas at 
ages younger than 2 Gyr. Further investigations are definitely re- 
quired to explain the metallicity gradients reported in some of 
the classical dSphs. 
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Fig. 16. Evolution of the metallicity with time. Top, model with a total mass of 3 x 10^ Mq. Bottom, the one with 9.5 x 10** Mq. The 
maps corresponds to a thin slice of 0.1 kpc. The colors code [Fe/H], from -4 to 0. 



Appendix A: Metals ejection and energy 
conservation 

In this appendix, we will describe precisely how the velocities 
of the neighboring stellar particles / are modified in order to im- 
prove the conservation of energy . 

The mass received by the particle j is : 

m'j - nij + Wij Me. (A.l) 

where is the mass ejected by the stellar particle /. Similarly, 
the mass of element k is : 



m'j(k) = mj(k) + w/j Ms(k). 



Here w,y is : 



Wii = 



mjWinj, hi) 



Pi 



The final mass the the stellar particle / is : 

m'j — nil ~ Mg. 



(A.3) 



(A.4) 



As the mass of particles changes, in order to conserve the 
(A. 2) total energy, it is necessary to modify the velocities of the parti- 
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cles. We neglect to correct the change in potential energy. Before 
the mass redistribution, the energy of the particles involved is : 



Zl 2 1 2 

^.v,. + -m,.v,., 



(A.5) 



and after, it is : 



where we have assumed that the particle / does not change its 
velocity. If v'J are chosen such that : 



1 ,2 1 



,2 1 2 1 



(A.7) 



the the summing over j gives : 



which implies that : 
E' = E. 



(A.9) 



The new velocity for each particle is deduced from Eq. |A.7 

,2 mjv] + WijM^v] rrij ^ Wij ^ 



' nij + Wij Me 



m. 
J 



(A.IO) 



The modification of the square of the velocity is due to two 
terms. The velocity is first decreased, due to the increase of the 
mass of the particle. The second term corresponds to the de- 
crease of energy of particle / due to its decrease in mass. As in 
practice, Wjj Mg is much smaller that nij, the change in veloc- 
ity is small. Only the norm of the velocity is affected during the 
modification of the velocities. 

The effect of this correction on the total energy is estimated 
by running the same simulation used in Fig. [3j but with the cor- 
rection switched off. The comparison of the evolution of the rel- 
ative energy is given in Fig. A. 1 The improvement is of about 
20% over 7 Gyrs, while its increase in CPU time is insignificant, 
we used it. 

We didn't observed any improvement in the linear momen- 
tum conservation. This is related to the poor conservation of the 
treecode method that does not fulfill Newton's third law and 
dominates the error 
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Fig. B.7. Effect of the star formation parameter c^,. The parameters for the initial conditions are Mtot = 7 x 10^ Mq, Pc.gas 
0.053 mn/cm^, rmax = 8 kpc. The parameters for the supernova feedback is esN = 0.05. 
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Fig.B.8. Effect of the star formation parameter psfr. The parameters for the initial conditions are Mtot = 3.5 x 10^ Mq, pc 
0.025 mn/cm^, rmax = 8 kpc. The parameters for the star formation and supernova feedback are c* - 0.05 and esN - 0.05. 
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Fig. B.9. Effect of the star formation parameter psfr. The parameters for the initial conditions are Mtot - 1 x 10^ Mg, 
0.053 mn/cm^, rmax = 8 kpc. The parameters for the star formation and supernova feedback are c* - 0.05 and esN - 0.05. 
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Fig. B. 10. Effect of varying the number of stars formed per gas particle A^*^. The parameters for the initial conditions are Mtot = 
5x10^ Mq, Pc „as = 0.037 mn/cm^, r^iax = 8 kpc. The parameters for the star formation and supernova feedback are c* = 0.025 and 
esN = 0.02. 
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Fig. B. 11. Effect of varying the IMF. The parameters for the initial conditions are Mtot = 2 x 10^ Mq, pc 
'"max = 8kpc. The parameters for the star formation and supernova feedback are c^, = 0.05 and esN - 0.03. 
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Fig.B.12. Effect of varying the IMF. The parameters for the initial conditions are Mtot = 8 X 10^ Mq, Pcgas 
'"max = 8kpc. The parameters for the star formation and supernova feedback are c* = 0.05 and esN - 0.03. 
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Fig. B. 13. Effect of varying the number of neighbors, A^ngb- The parameters for the initial conditions are Mtot = 4.5 x 10** Mg 
Pc.gas — 0.033 mn/cm^, rj^ax — 8kpc. The parameters for the star formation and supernova feedback are — 0.05 and egN — 0.01. 
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Fig. B. 15. Effect of the total mass. We vary r^ax, keeping a constant central gas density Pcgas = 0.044 niH/cm^. The parameters for 
the star formation and supernova feedback are =0.1 and esN - 0.05. 
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Fig. B.16. Effect of varying the central gas density, for a constant total mass Mtot = 1x10** Mq. The parameters for the star formation 
and supernova feedback are =0.1 and esN = 0.05. 
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Fig. B.17. Effect of varying the central gas density, for a constant total mass Mtot = 5x10** Mq. The parameters for the star formation 
and supernova feedback are =0.1 and esN = 0.05. 



